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ABSTRACT 

We perform a weak-lensing study of the nearby cool-core galaxy clusters, Hydra A 
{z = 0.0538) and A478 (z = 0.0881), of which brightest cluster galaxies (BCGs) host 
powerful activities of active galactic nuclei (AGNs). For each cluster, the observed 
tangential shear profile is well described either by a single Navarro-Frenk-White model 
or a two-component model including the BCG as an unresolved point mass. For A478, 
we determine the BCG and its host-halo masses from a joint fit to weak-lensing and 
stellar photometry measurements. We find that the choice of initial mass functions 
(IMFs) can introduce a factor of two uncertainty in the BGG mass, whereas the BGG 
host halo mass is well constrained by data. We perform a joint analysis of weak- 
lensing and stellar kinematics data available for the Hydra A cluster, which allows 
us to constrain the central mass profile without assuming specific IMFs. We find 
that the central mass profile (r < 300 kpc) determined from the joint analysis is in 
excellent agreement with those from independent measurements, including dynamical 
masses estimated from the cold gas disk component. X-ray hydrostatic total mass 
estimates, and the central stellar mass estimated based on the Salpeter IMF. The 
observed dark-matter fraction around the BGG for Hydra A is found to be smaller 
than those predicted by adiabatic contraction models, suggesting the importance of 
other physical processes, such as the the AGN feedback and/or dissipationless mergers. 

Key words: galaxies: clusters: individual (Hydra A, A478) - gravitational leasing: 
weak 


* Based on data collected at Subaru Telescope, which is operated 
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1 INTRODUCTION 


The standard cold-dark-matter (CDM) paradigm excel¬ 
lently describes the statistical properties and growth 
of cosmic structure from the early universe to the 
present day epoch, as observed fr om the cosmic mi¬ 
crowave background radiation ( e.g. iHinshaw et alJ l2013l : 
IPlanck Collaboration et alj 20lj) a nd lar g e scale cluster¬ 
ing o f galaxies (e.g. iPercival et ^ I 2001 I : lAnderson et alJ 
|20141 . Cosmological simulations of collisionless CDM pre¬ 
dict that the quasi-equilibrium mass density profile of dark- 
matter halos formed by collisionless gravitational dynam¬ 
ics is n early universal and self - simila r over wide mass 
ranges dNavarro. Frenk fc Whitd Il996l . hereafter NFW). 
The spherically-averaged density profile exhibits a soft cusp 
{p oc in the center and has steeper slopes at larger 

halo radii (p oc in the outskirts). Recent observations 
of stacked weak gravitatio nal lensing signals for massive 
galaxy clusters at z ^ 0.2 dOgu ri et al.| 201^Okabee£_^ 
2 OIOI. 2013l:l0kabe fc Smithll2015l:IUmetsu et al.ll201ll.l20l4 


2015b nNiikura et al.ll2015l f revea ed that the observed clus¬ 

ter lensing profiles are in excellent agreements with the NFW 
density profile predicted for collisionless CDM halos. How¬ 
ever, since the innermost cluster radius for weak lensing is 
typically limited to r 100h“^ kpc due to the finite reso¬ 
lution of weak-lensing observations, it is practically difficult 
to directly constrain the innermost cluster mass distribution 
around bright cluster galaxies (BCGs). 

The total mass density profiles of clusters in the cen¬ 
tral region have been a subject of intense theoretical and 
observational studies. Numerical simulations of collision¬ 
less dark matter show that the central c usp slope of dark - 
matter halos is not strictly self similar dCao et al.l I 2 OI 2 II . 
Baryonic effects cannot be ignored in the central region 
where the cooling time is shorter than the Hubble time 
and rich stellar populations within BCGs are abundant. 
Stars form from radiative gas cooling, and subsequently 
dark matter is pulled inward in response to the central 
condensation of baryons. A cuspy dark-matter density pro¬ 
file is formed through such energy exchange processes be¬ 
tween dark matter and b aryons, which is refereed to as 


_ ( e.g. 

iBlumenth^ et al l ll98d : lGao et al.ll2004l : lGnedin 3t al 


adiabatic contraction (e.g. |EggemHj^ndemBe]l_&_Sandnge 


2004 1. On the other hand, if BGGs are formed by the accre¬ 
tion of clumpy structures such as satellite galaxies and glob¬ 
ular clusters, a cored density profile is expected to form by 
dynamical friction (e.g.|El-Zant. Shlosman fc Hoffmarj200ll : 

iLackner fc Ostrikerj I 2 OI 0 I I. Gas outflows driven by active 
galactic nuclei (AGNs) can suppress condensation of gas and 
stars, resulting in flatten ing of the dark-matter distribution 
in th e central region (e.g.|Teygaer_et__ali201l|; Martizzi et al.l 


I 2 OI 2 I : iRagone- Figueroa" Granato fc Abadill2012ll . 

Stellar kinematics of BGGs is one of the most power¬ 
ful means of measuring t he total mass profile in the cen- 
tral region of galaxies (e.g.lGappellari et al.ll200^ f. Recently, 
I Newman et al'l (l2009l . 1201^ 2013li measured the total mass 
profiles of relaxed clusters and decomposed them into into 
stellar and dark-matter components from joint measure¬ 
ments of strong/weak lensing and stellar kinematics within 
BGGs, under the assumptions of dynamical equilibrium, 
isotropic orbits, and const ant stellar mass-to-light ratio. 
lOguri. Rusu fc Falcd (l2014l l carried out a joint ensemble 


analysis of strong lensing and stellar photometry data for a 
large sample of 161 elliptical galaxies to measure their aver¬ 
age mass density profile. Stellar mass estimates for galaxies 
can be obtained from stellar photometry using scaling rela¬ 
tions between galaxy luminosity and stellar mass, whereas 
the stellar initial mass function (IMF) governs the normal¬ 
ization of the stellar mass to luminosity relation. Hence, the 
uncertainty of the stellar IMF can give rise to a systematic 
bias in stellar mass estimates. In order to break the degener¬ 
acy between the IMF uncertainty and the inner s lope of the 
mass density profile, lOguri. Rusu fc Falcd (120141') proposed 
to use quasar microlensing observations. 

We have two ide al clusters. Hydra A (2 = 

0.0538; ISato et al.l l2012l'l and Abell 478 (2 = 0.0881; 


iMochizuki et al] ~ 2014l ~).~ to study the mass distributions 
within both clusters and BGGs. The two clusters have 
been selected and observed as tar gets of joint Subaru weak- 
lensing and Suzaku X-ray studies ( Okabe et al .l2014bl) . They 
are well-known cool-core clusters hosting central AGNs 
insid e their BGGs. Recent X-ray and radio o bservations 
(e.g. j McNjunaxaeL^ 200^David_et_ah 200 ll; Wise et all 


20T^ Sun et al.l2003l : ISanderson. Finoguenov fc Mohj|2005l : 


Diehl et al.l 2008l l have shown that giant radio lobes trig¬ 
gered by prominent AGN jets are interacting with the sur¬ 
rounding intracluster medium (ICM). Isobaric cooling time 
scales in the central region are shorter than 1 Gyr. Hence, the 
clusters provide us with unique environments for studying 
the effects of baryons acting on the dark-matter distribution 
in the cluster center. 

Since these clusters are at very low redshifts, their large 
angular extents enable us to detect weak-lensing signals 
around the BCGs, in particular for Hydra A, and to measure 
the radial mass distribution, by using the unique combina¬ 
tion of weak lensing and stellar kinematics or stellar pho¬ 
tometry. For these very nearby clusters, especially for the 
Hydra A cluster, strong-lensing information is not needed 
to resolve the central matter distribution. 

This paper is organized as follows. In Section [5] we 
conduct shape measurements, examine the redshift depen¬ 
dence of the contamination level by member galaxies, and 
securely select background galaxies. We measure cluster 
masses in Section [S] perform joint analyses in Section [I] 
and discuss the results in Section [5] Section [6] is devoted 
to the summary. A comparison of weak-lensing masses and 
Suzaku X-ray obser vables is given in our previous paper 
dOkabe et al.ilioidlj) . We use Ilm.o = 0.27, Da = 0.73 and 
Ho = 70/170 kms~^ Mpc“^. 


2 DATA ANALYSIS 

2.1 Data Analysis and Shape Measurement 

We conducted observations of Hydra A and Abe ll 478 using 
the Subaru/Suprime-cam llMivazaki et al.ll2002lf . in i' and 
g' bands, on Jan 7th, 2013. The i imaging was taken to 
measure ellipticities of galaxies for the weak-lensing shape 
analysis. The two-band photometry in the g' and i' bands is 
to minimize contamination of the background source cata¬ 
log by unlensed cluster member galaxies. The Hydra A and 
Abell 478 fields are covered by a mosaic of 8 and 9 point¬ 
ings, respectively, as shown in Figure [T] We note that, for 
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Table 1. Summary of Subaru observations for Hydra A and Abell 
478. Each pointing is shown in Figure [T] The seeing FWHM is 
for the i'-band used for shape measurements. 


Name 

i> 

9' 

Seeing 


[min] 

[min] 

[arcsec] 

HYDRAA.00 

15.0 

11.7 

0.75 

HYDRAA.02 

15.0 

10.0 

0.71 

HYDRAA.10 

15.0 

11.7 

0.73 

HYDRAA.il 

15.0 

11.2 

0.71 

HYDRAA.12 

15.0 

10.0 

0.63 

HYDRAA.20 

15.0 

11.7 

0.67 

HYDRAA.21 

15.0 

11.2 

0.71 

HYDRAA.22 

15.0 

10.0 

0.65 

ABELL478_00 

15.0 

11.7 

0.67 

ABELL478.01 

15.0 

11.7 

0.79 

ABELL478.02 

15.0 

11.7 

0.73 

ABELL478.10 

15.0 

11.7 

0.65 

ABELL478.il 

15.0 

11.7 

0.79 

ABELL478.12 

15.0 

11.7 

0.71 

ABELL478.20 

15.0 

14.0 

0.69 

ABELL478.21 

15.0 

11.7 

0.69 

ABELL478.22 

15.0 

11.7 

0.71 


correlations is consistent with zero. The reduced shear sig¬ 
nal is obtained by applying a correction for the isotropic 
smearing effect as ga = {Pg^)ape'p- Here, e„ is the image 
ellipticity after the anisotropic PSF correction and is 
the pre-seeing shear polarizability tensor. Since the measure¬ 
ment of P^p is very noisy for faint individual galaxies, P®^ is 
calibrated using galax ies detected with hig h signal-to-noise 
ratio > 30, following lOkabe et ah (|2014a ). Sim i lar cal ibra¬ 
tion procedures ar e developed by lUmetsu et al.l (l2010ll and 
lOguri et al.l (l2012h. Our shear calib ration on Qa results in 
2-3% accuracy ( Okabe et al.|[^14ah . 

We measure magnit udes of galaxies using SExtractor 
llBertin fc Arnout^ Il996l l. The SExtractor configuration is 
optimized for shape measurements of faint galaxies. For 
each object, we compute a total Kron-like magnitude and 
an aperture magnitude for color measurements in the AB- 
magnitude system. To perform color measurements, the 
f'-band data are degraded to the PSF of g'-band data. 
For aperture photometry, the aperture diameter is set to 
1.5 times the worst-seeing FWHM. Finally, we combine 
SExtractor-based photometry with weak-lensing shape mea¬ 
surements to create a weak-lensing-matched photometry 
catalog. 


2.2 Background Selection 


Hydra A, the middle-east pointing was not available as lim¬ 
ited by the number of allocated nights and bright stars in 
that region. The observing conditions are summarized in 
Table m The typical seeing is ~ 0.7” in the i' band for 
weak-lensing shape measurements. The exposure times for 
the i' and g' bands are ~ 15 and ~ 12 min, respectively. 
These are about half of those typically used in our previous 
weak-lens ing studies of intermediate-redshift clusters at z ~ 
0.2 (e.g., Ofabe k. Umetsul l2008l : lOkabe et al.l l2010l . l2013l : 
lOkabe fc SmithI 2015lj . We use d the standard Suprime-Cam 
reduction software SDFRED dYagi et al.ll2002l : lOuchi et al.l 
|2004^ modihed to accommodate new CCD chips, for flat- 
helding, instrumental distortion correction, differential re¬ 
fraction, point-spread-function (PSF) matching, sky sub¬ 
traction and stacking. An astrometr ic calibration was con - 
ducted with 2MASS point sources (ISkrutskie et al.l [200611 . 
Typical residual astrometric offsets are no larger than the 
CCD pixel size. 

We determine the mass distribution in nearby clus¬ 
ters using weak gravitationa l leasing. Weak lensing anal¬ 
ysis i s performed following lKalserj_S 2 uires_&^roadhur^ 
1199^ with sorn e modihcations l Okabe et al.ll2013l . 12014'^ 
Okabe_^_Sinlth|_[201^ . Technical details are described in 


Okabe et al.l (l2014al l. The image ellipticity Ca of an object 
detected in the i' band data is measured from weighted 
quadrupole moments of the surface brightness distribu¬ 
tion. The PSF anisotropy is corrected for using second- 
order bipolynomial functions of the stellar anisotropy ker¬ 
nel, which is expressed with KSB’s smear polarizability ten¬ 
sor and image ellipticity of unsaturated stars. To examine 
the validity of anisotropic PSF corrections, we have checked 
residual syste matics in corrected e llipticities as described 
in Appendix of lOkabe et al.l (l2014al') . Two-point correlation 
functions between star and galaxy ellipticities after the cor¬ 
rection are of the order of signihcantly improved 

from before the correction. This level of residual 


It is of prime importance for an accurate mass measure¬ 
ment to securely select background galaxies. Contamina¬ 
tion of unlensed member and/or foreground galaxies in 
the shear catalog leads to a systematic underestimation 
of cl uster lensing mass, ref e rred to as a dilution effect 


Broadhurst et al. 20051: Umetsu fc BroadhurstI 2008|; 


Okabe et al.]|20ld. 1201 iL l2013l . l2014all . 


The dilution effect is increasing as a projected cluster-centric 
radius is decreasing, because the fraction of member galax¬ 
ies to background galaxies is increasing. The effect is more 
dominant for more massive clusters. We select background 
galaxies in color-magnitude plane, f ollowing studies for clus¬ 
ters at 2 ^ 0.2 or highe r redshift (iBroadhurst et aT1l2005l : 
lOkabe et al.ll2olol . l2013l l. A majority of unlensed galaxies 
to dilute lensing signals are faint and small member galax¬ 
ies around the red sequence and apparent magnitudes i' ~ 
22 — 26 ABmag in the color-magnitude plane. Although the 
red sequence of bright member galaxies are clearly identi- 
hed in color-magnitude plane because of a passive evolution, 
the color distribution for faint galaxies becomes unclear as 
the apparent magnitude becomes fainter. It is therefore very 
difficult to discriminate between member and background 
galaxies in the faint end. We thus investigate the color dis¬ 
tribution of member galaxies using lensing signals in order to 
minimize the contamination. We measure a mean tangential 
distortion strength in the central region of clusters by chang¬ 
ing lower limits of color for the background selection, {gO{> 
AC) = Y.i9+A> ^C)wg,i{> AC)/Y..Wgp{> AC). Here, 
AC = (p' - i') - {g' - i') E/go is as a function of color offset 
from the red sequence, g+p is the tangential component of 
reduced shear of the i-th galaxies with respect to the cluster 
center, given by 

g+,i = -(pi,iCos2(p-|-p2.isin2(/3), 

and Wgp = XjtPgp + o?) is a statistical weight. Here, ip 
is the position angle between the Erst coordinate axis on 
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Figure 1. Pointing maps for Hydra A (left) and Abell 478 (right). The horizontal and vertical axes are R.A. and Decl. offsets from BCGs, 
in units of arcmin. For each cluster, the circle represents the position of the BCG: (139.524°,—12.095°) for Hydra A and (63.355° ,10.465°) 
for A478. The pointing IDs ('TableFTTl are shown for each cluster. The adjacent pointings have an overlap of 2 arcmin. The hatched region 
indicates that there is no data for the present study. 


the sky and the vector connecting the cluster center and 
the galaxy position. The position of the BCGs is adopted 


timate and a is the softening constant variance represent¬ 
ing t he scatter due to the intrinsic ellipticity of the galaxies 


mg t ne scatter Que to tiie intrinsic ellipticity oi the galaxies 
(e.g., Hoekstra. Franx Sz Kuiikenlj200o f:lHarnan a et al.ll 20Qd : 


lOkabe et alj 120101 : lUmetsu et al]l20ldnOguri et al.ll2012ll . 

We choose a = 0.4. 


In order to clarify the dilution effect in noisy lensing sig¬ 
nals, the mean distortion strength is cumulatively computed 
as a function of color offset from the red sequence, AC. We 
choose the central region of 100 — SOO/i^q^ kpc from brightest 
cluster galaxies (BCGs). We also compute the mean lensing 
depth (Dis/Ds). Since it is difficult to estimate photometric 
redshifts of individual objects using two bands, we used the 
COSMOS photometric redshift catalog (lllbert et al.ll200'^ l 
estimated by combining 30 bands. The mean photomet¬ 
ric redshift is computed with a statistical weight of Wg by 
matching with the COSMOS photometric redshift in color- 
magnitude plane. If there were no contamination of unlensed 
member galaxies, it is expected that a “Dilution” estimator, 
which is defined by the ratio between the mean lensing sig¬ 
nal and depth, D = {g+)/ {Dis/Da){> AC), is constant irre¬ 
spective of the lower limit of the color offset, AC. Here, Dis 
and Ds are the angular diameter distance between cluster 
(lens) and source redshifts and between observers and source 
redshifts, respectively. Therefore, D is a good estimator to 
investigate the color distribution of member galaxies. 

The resultant D is shown in the top panel of Figure [21 
The red diamonds and blue circles denote D for A478 and 
Hydra A, respectively. Although the error is large, D of A478 
is constant at AC ~ 0, increases with AC at AC ~ 0.1, and 
becomes flat at AC > 0.3. The small D at AC ~ 0 indi¬ 
cates a presence of member galaxies. On the other hand, D 


of Hydra-A shows a random distribution for AC, suggest¬ 
ing no significant feature of the contamination. We model a 
color distribution of D, given by A+B tanh((AC—ACo)/^). 
Here, A is the normalization at AC = 0, A + B is a mean 
lensing strength of pure background galaxies, and the sec¬ 
ond term represents a constant D at AC 0 and AC > 0.3. 
As shown in Figure (21 the best-fit model (red solid line) of 
A478 well describes the data. As for Hydra A, the best- 
fit model is consistent with no contamination. We derive a 
color distribution of contamination level of unlensed galax¬ 
ies in background shear catalog (the bottom panel of Figure 
0 . The maximum contamination level is at most 17% at 
AC = 0. In order to investigate a cluster-redshift depen¬ 
dence of the dilution effect, we compute the mean D within 
the same physical radius for 50 clusters at z ^ 0.2 from the 
Local Cluster Substructure Survey (LoCuSS). Here, AC is 
the color offset from the red sequence in the i! — V band. The 
ca lculation and mode l for stacked 50 clusters are described 
in lOkabe et al.l (l2013l L The maximum contamination level 
for 50 clusters (green dashed line in the bottom panel) is 
~ 40% at AC = 0. The virial mass of A478 is more mas¬ 
sive than the average of 50 clusters (Section[2| and our data 
is shallower than those of 50 clusters, but nonetheless the 
contamination is less than half of those of clusters a ~ 0.2. 

There are two reasons for low contamination of the 
very nearby cluster {z ~ 0.09). First, since the colors of 
red-sequence galaxies become more blue as the redshift de¬ 
creases, the number of red background galaxies increases. 
Second, the ratio of member to background galaxies for very 
nearby clusters becomes lower because an apparent covering 
area becomes larger. Therefore, the dilution effect for clus¬ 
ters at lower redshifts becomes drastically less significant. It 
is consistent with the feature of the very nearby cluster of 
Hydra-A {z ~ 0.05) . 
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Figure 2. Top: dilution estimator, D = {g+)/ {Dig / Ds), as a 
function of color offset from the cluster red sequence. The red di¬ 
amonds, blue circles, and green crosses show D for A478, Hydra A, 
and an ensemble of 50 clusters from the LoCuSS project, respec¬ 
tively. The relative lensing strengths D for A478 and the LoCuSS 
50 clusters increase with increasing AC, while D for Hydra A is 
randomly distributed. The solid lines represent the best-fit mod¬ 
els to describe the color distribution of member galaxies. Bottom: 
contamination level calculated by the respective best-fit model as 
a function color offset. The red-solid and green-dashed lines show 
the results of A478 and the LoCuSS 50 clusters, respectively. 


In order to securely select backgro und galaxies, we em - 
ploy 1% contamination level following lOkabe et al.l (|201.‘ll l 
and define background galaxies with AC > 0.17 for A478. 
As for Hydra A, we do not apply the color cut because of 
the tiny contamination level, because D is almost constant in 
the wide color range of — 1 < AC < 1. In order to assess our 
background selection for Hydra A, we make another back¬ 
ground catalog selected with the same color cut of A478, 
conduct a weak lensing mass measurement (Section |3]) and 
confirm the result dose not change. The background number 
densities for A478 and Hydra A are ribkg = 6.2 [arcmin“^] 
and 28.9 [arcmin“^], respectively. The mean source redshifts 
for A478 and Hydra A are (zs) — 0.643 and (zs) — 0.591, 
respectively. 


3 MASS MEASUREMENTS 

Tangential distortion signals contain full information on the 
lensing signals of clusters. Since clusters can be considered 
as gravitationally bound spherical objects at zeroth order, 
a measurement of tangential shear profile is a powerful way 
to measure cluster masses. The tangential shear is measur¬ 
able by averaging over sufficient background galaxies, in or¬ 
der to suppress shape noise attributable to random intrin¬ 
sic ellipticities of background galaxies. We measure a mean 
tangential shear (p+.i) as an average of the tangential shear 
component of background galaxies residing in the i-th radial 


bin with the statistical weight of Wg. The other distortion 
component, gx, which is the 45 degree rotated component, 
is also estimated. Hereafter, we refer to {g+(x),i) as g+(x),i- 
In order to interpret the lensing signals, we here briefly 
introduce mass models as the cluster mass density profile. 
The reduced tangential shear prohle for mass models can be 
expressed by 


9+0 


k{< 6) — K{d) 
1 - k(6») 


( 1 ) 


where R and k are the mean surface mass density within 
a radius 9 and the surface mass density at 9, respectively. 
The lensing signals are computed by R and k for given mass 

models. __ 

The NFW model llNavarro. Frenk fc Whitelll996l . Il997l 'l 
is a theoretically well-motivated mass model according to 
numerical simulations of dark matter particles. They found 
that the matter density profile of CDM halos is well de¬ 
scribed by an analytic function 

PNFw(r) = -i—f —IFTT— 1 —11 ’ O 

{r/rs)(l + r/rsY 

where ps is the central density parameter and is the scale 
radius. The asymptotic mass density slopes are p oc r~^ and 
r~^ for r <^rs and r ^ Ts, respectively. The NFW profile is 
specihed by two parameters of Ma which is a mass enclosed 
within a sphere of radius ta, and the halo concentration 
CA = TA/rs. Here, ta is the radius inside of which the mean 
density is A times the critical mass density, pcr{z), at the 
redshift, 2 . For a comparison, we model a singular isothermal 
sphere (SIS) profile for the main halo, given by 


Psis(r) 


2nGr^ ’ 


(3) 


where crj is a one-dimensional velocity dispersion. We also 
consider the BCG mass model to estimate the lensing signals 
for the BCGs. Since the innermost radii for the profiles are 
outside the BCGs, the tangential shear from the BCGs is 
sensitive to the interior mass within 9, in other words, n{9) = 
0 in Equation O- Since the internal structure for the BCGs 
cannot be resolved by lensing signals, we reasonably employ 
the point mass, Mpt, to estimate the BCG mass, in terms 
of 9+0 oc Mpt/(7r61^). 

Given specific mass models, we perform fit fo the 
tangential shear profile, g+. The is defined as followed: 


2 

Ag 


E ( model\ ( model\ 

\9+,i - 9+,i ) Cig {g+,j - g+^g ) , 


(4) 


where 


Cij — Cstat,ij C+SS,i, 


(5) 


The first term of C is the statistical noise relevant to the 
intrinsic ellipticity of background galaxies. The statistical 
noise for the shear is given by Cstat.ij = Pg+Sijl'^j where 
(Tg,i is the statistical uncertainty of the two components in 
the i-th radial bin and 5ig is Kronecker’s delta. The sec¬ 
ond term, C+sB+j, is the error covariance matrix of uncorre- 
lated large-scale struc t ure (LSS) along the line-of-sight, (e.g . 
Schne ideiieUalJ ]^98|; Huekstra 200^ Umetsuet_aL 2011, 


20 14 Oguri et al.ll201C : IOguri fc Takadall201ll : lOkabe et ahl 
20131 . 2014al l. estimated from th e weak-lensing power spec¬ 
trum with WMAP7 cosmology (iKomatsu et al.fl201ll h The 
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cosmic shear correlation function leads to non-zero off- 
diagonal elements for Cnss.ij- 

The signal-to-noise ratio (S/N) for the tangential shear 
profile is computed by 

(f) =Yl9+,iC^r/g+J- (6) 

ij 

3.1 A478 

The tangential reduced-shear profile for A478 is shown in 
Figure [3] in a range of lOOkpc ^ r ^ 5.5Mpc. The inner¬ 
most radius, determined by requiring a sufficient number of 
background galaxies, is larger than the effective radius of 
the BCG. The S/N of the tangential shear profile is 8.1. 
The tangential shear profile exhibits a clear curvature. On 
the other hand, the B-mode, gx, shows both negative and 
positive values, which is consistent with null signal. We plot 
6 ■ gx (0) in the bottom panel of Figure [4] so that the scatter 
is independent of rad ius for logarithmically spaced binning 
llMivatake et al.ll2013l l. 

We fit a single NFW model or two components of the 
NFW model as the smooth matter density profile and the 
point source {g+ = (?nfw + ffpt ; NFW-fpoint) to the tan¬ 
gential shear profile (Table [2]). The minimum chi-square 
are Xmin(d.o.f) = 2.18(7) and 2.12(6) for the first and 
second model, respectively, where the degree of freedom 
is presented in the round bracket. Here, since errors are 
mainly attributable to the variance of intrinsic ellipticity, 
~ 0.4, in addition to measurement errors, Xmin be¬ 
comes small. It shows that two models are acceptable as 
the cluster mass model. In order to estimate the validity of 
the additional point source, we calculated F test probabil¬ 
ity, 0.7, indicating that there is no strong reason to add the 
point source to describe the current data. The best-fit SIS 
profile ((T„ = 958.67lg2'ge kms“^) is deviated from the cur¬ 
vature of the tangential shear profile , albeit still acceptable 
(xLn(d.o.f) = 9.13(8)). 

3.2 Hydra A 

A large apparent size of Hydra A allows us to estimate 
the tangential shear profile in a wide radial range of 
30 kpc ^ r ^ 3 Mpc, as shown in Figure (4] In particular, 
lensing signals at r 50 — 100 kpc is detected solely by 
weak-lensing analysis. It is an advantage of weak-lensing 
study of very nearby clusters, because cluster central re¬ 
gions at higher redshifts are in the strong-lensing regime. 
For instance, an innermost radius of a stacke d lensing pro¬ 
file for 50 clusters z ~ 0.2 dOkabe et al.|[201^ is 70h/Q^kpc. 
The tangential profile is highly detected at S/N = 7.0. The 
profile shows two characteristic features : a sharp decrease 
in the central region (r < 100 kpc) and a clear curvature in 
the range of 100 kpc ^ r ^ 3 Mpc. The B-mode values, gx, 
are randomly distributed cross the null signal. 

First, we fit the profile with the NFW model and obtain 
a high concentration parameter Cvir = 9.37/l3 g 4 with virial 
mass Mvir = 2.74l4’2i 10^“^ Mq. The best-fit mass profile 

shows a clear curvature (Figure [4|). However, we found that 
the best-fit model is lower than the observed lensing signals 
at r 50 kpc and r ^2 Mpc. It suggests that the virial mass 
is underestimated because the mass estimate is sensitive to 





1 10 50 

9 [arcmin] 


Figure 3. Radial profiles of the tangential reduced-shear com¬ 
ponent (top panel), and the 45°-rotated component (bottom 
panel), d-gx , for A478. The vertical error bar of each data point (f- 
th bin) shows the diagonal element of the error covariance matrix, 
(Equation [Sjl , for the top panel and for the bottom 

panel, respectively. The tangential reduced-shear profile shows a 
clear radial curvature, in agreement with the steepening NFW 
profile. The red-solid, green-dotted, and blue-dashed lines show 
the best-fit NFW, NFW-|-point, and SIS profiles, respectively. 


lensing signals in the outskirts. A similar deviation is found 
in the best-fit SIS model = 631.93'^*Y'’rr kms“^). On the 
other hand, when we fit the NFW model to the lensing pro¬ 
file beyond 150 kpc, the best-fit concentration becomes lower 
Cvir = 4.52l4 ®g (Table [2]). The best-fit profile well describes 
the data at r ^ 300 kpc, but significantly underestimates 
the signals at r 100 kpc. 

It indicates that the inner slope of the tangential profile 
becomes stepper than expected from the outskirts. The steep 
inner slope suggests one possibility that an extra component 
exists inside the innermost radius. In particular, the slope 
for the lensing profile at a few inner bins is proportional 
to ~ —2. We thus consider two components of the NFW 
model and the point mass (NFW-l-point). The best-fit pa¬ 
rameters are shown in Tabled The Xmin(d-o.f) is improved 
from 8.62(18) to 4.02(17). The best-fit profile well describes 
the data in the full radial range (Figures m. In order to 
estimate the validity of the additional point source, we con¬ 
ducted F test and found that the probability 4 x 10~® is sig¬ 
nificantly small. It is thus reasonable to add the extra com¬ 
ponent to the smooth mass component for describing the ob¬ 
served lensing signals. In other words, the two-components 
model better describes the tangential profile. 
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Table 2. Mass estimates for Hydra A and A478. The NFW mass {M^^) and the point mass (Mpt) are in units of 10^“^ Mq and 
10^^ respectively. 


Cluster 

model method 

mNfw 

vir 

mNfw 

mNfw 

^NFW 

^^^1000 

/vtNFW 

-^‘^2500 

-^pt 

A478 

NFW 

16.381®-™ 

12.931®;^® 

8-67l?;i^ 

5 . 951 ?-?? 

0 1-7+0.72 

0.1 / _o.72 

— 

A478 

NFW-l-point 



8.561?-™ 

5.74l?-®5®3 

2 93+1-28 

< 25.90 

Hydra A 

NFW (r > 150 kpc) 

4.661?;™ 

31 

2 

^•^U_0.76 


U.»0_Q 28 

- 

Hydra A 

NFW+point (full radial range) 

4.88l®;29 

3.661?-®? 

Q Qy-|-0.98 

-0.77 

1 43+0-^^ 

0 64'*”^'^^ 

z. -,r»+2.47 

6-1012.50 


r [hyo^Mpc] 



1 10 50 


9 [arcmin] 


Figure 4. Radial profiles of the tangential reduced-shear com¬ 
ponent (top panel), and the 45°-rotated component (bottom 
panel), 6 ■ gx, for the Hydra A cluster, shown over a wide radial 
range of 30kpc r 3Mpc. The red-solid, magenta-dot-dashed, 
green-dotted, and blue-dashed lines show the best-fit NFW pro¬ 
files derived from the data at r > 30kpc and r > 150 kpc, the 
best-fit NFW-l-point and SIS profiles, respectively. 

4 JOINT ANALYSES 

4.1 Joint Weak-lensing and Stellar Photometry 
Stndy of A478 

We study the mass distribution of A478 using weak-lensing 
and BCG (PGC014685) photometry measurements. 

The tangential shear profile (Section 13.111 for A478 is 
well described by a single NFW model or a two-component 
mass model including the BCG as an unresolved point mass. 
Gomplementary information at small scales is useful to di¬ 
rectly determine the BGG mass profile as well as to better 
constrain the total mass profile in combination with lensing 
data. 

For this purpose, we use stellar mass estimates of the 
BGG from a photometry. We estimate the stellar mass for 
the BGG using A-band photometry using the A-correction 
and galactic extinction from 2MASS extended source catalog 
llSkrutskie et al.ll2006l ). Since the A-band luminosity is sen¬ 
sitive to the light of old stars and less sensitive to recent star- 


formation activities, it serves as a reasonable proxy for stel¬ 
lar mass. The difference in A-band transmission between the 
2MASS and tra ditional Johnson A filters is not important 
for our purpose (| Garpen terj|200lh . We e mploy both empirical 
( [Ajnouts et al.l [20071 ) and theoretical llLonghetti fc Saraccd 
I 2 OO 9 I) scaling relations between the stella r mass and A - 
band lum inosity. We con sider the SalpeterjjSalpetej l 19551 ). 
Chabrier llGhabrierll2003l) and Kroupa (lKroupal|200ll) IMFs. 
The mass uncertainty is estimated by taking into account 
the photometric erro rs, 15% intrinsic scat ter for the empir¬ 
ical scaling relation (lArnouts et al.l l2007l ) , and systematic 
uncertainties from different scaling relations. 

We assume a pseudo isothermal mass distribution 
(PIEMD) mode l (lElia sdottir et al.l l2007l : iLimousin et al.l 
I2OO9I : iNataraian et al.ll2009l) to describe the stellar mass dis¬ 
tribution of the BCG. The PIEMD model is specified with 
three parameters, namely, the core radius (a), the scale ra¬ 
dius (s), and the normalization (po), in the following form: 

(1 + rya^Xl + rys^) ^ 

The projected sufrace luminosity profile for the BCG is 
shown in Figure (5] The profile is fitted with the PIEMD 
model convoluted with the PSF kernel. We find that the 
outer edge of the luminosity profile is slightly truncated and 
deviated from expectations for the PIEMD model. Although 
the truncated PIEMD model would be better to fit the over¬ 
all profile, the deviation is small for this study and thus we 
adopt the PIEMD model here. We obtain a = 2.7kpc and 
s = ll.Skpc. Assuming a scaling rel ation between cusp y ra¬ 
dius and the luminosity for BCGs l|Lauer et al.l 12004 ). the 
mean core radius is expected to be a ~ O.Skpc. The best- 
fit core radius is larger than the mean one calculated from 
the scaling relation, which might be caused by a limitation 
of ground-based telescope because the apparent size of the 
expected core radius, ~ 0.5arcsec, is small. We find the 
best-fit solution by minimizing the combined function, 

= Xg + As, with xt defined by 

2 (M* — Mpiemd,*)^ 

X. = ^-72-(8) 

"m. 

Here, M« and 5m, are the stellar mass and its uncertainty 
estimated with specific IMFs within the measurement radius 
17 kpc, respectively. Mpibmd is the model for the BCG. We 
determine the normalization of the PIEMD model and the 
NFW parameters from the joint fit. 

Figure |6] shows the best-fit model for the tangen¬ 
tial shear profile and the BCG stellar mass profile. The 
Chabrier and Kroupa IMFs give total stellar mass esti¬ 
mates of M§cg,* ~ ~ 8 X iC^Mq, compared to 
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Figure 5. Projected surface luminosity profiles for BCGs in A478 
(red circles) and the Hydra A cluster (blue squares) in arbitrary 
units. The solid and dotted lines denote the best-fit profiles for 
A478 and the Hydra A cluster, respectively 


M|cg,* = (15 ±2) X with the Salpeter IMF, where 

the errors mainly account for systematic uncertainties from 
different scaling relations. The IMF uncertainty gives rise to 
a s ystematic bias on stel l ar ma ss estimates, as pointed out 
by lOguri. Rusu fc Falcd (l2014h . Hence, in dependent con¬ 
straints, such a s from stellar kinem atics dNewman et^ 
200£ . l201ll . [2013l ) quasar microlensing dOguri. Rusu Sz Falco 
2014) observations, are needed to meaningfully break the 
IMF uncertainty (see Section 14.211 . Despite the IMF un¬ 
certainty, the best-fit NFW mass model (Mvir = 16.48 x 
10^^ Mq and Cvu = is not sensitive to the choice 

of IMF because the leasing signal from the main cluster halo 
is sufficiently high. We note that the result does not signifi¬ 
cantly change when we adopt a — O.Skpc. 


4.2 Joint Weak-lensing and Stellar Kinematics 
Study of Hydra A 

The tangential shear profile for the Hydra A cluster is well 
described by the two-component mass model of the BCG’s 
point mass and the smooth NFW profile (Section 13.211 . The 
smooth mass component is determined by fitting leasing sig¬ 
nals at outer cluster radii, whereas the point mass measure¬ 
ment is sensitive to lensing signals at innermost bins because 
the BCG lensing signal is em bedded in the lensing s ignals 
of th e main cluster. Recently iNewman et al.l ll2009l . l201ll . 
I2OI3I ) performed a joint analysis of weak/strong-lensing and 
stellar kinematics data to break the IMF uncertainty and 
measure the mass profile from the BCG to the cluster main 
halo. They demonstrated that the BCG stellar kinematics 
is powerful for determining the BCG mass profile and the 
cluster inner mass profile under two assumptions of dynamic 
equilibrium and mass tracing light. Following their pioneer¬ 
ing studies, we here conduct a joint weak-lensing and stellar- 


Figure 6. Tangential reduced-shear profile and BCG stellar mass 
profile (inset) for A478. The red-solid line represents the best- 
fit solution for the total lensing signal with the NFW (blue 
dot-dashed line) plus BCG (Salpeter IMF) model. The green- 
dotted and magenta-dashed lines show the best-fit lensing pro¬ 
files and BCG stellar mass profiles, obtained with the Salpeter 
and Chabrier IMFs, respectively. 


kinematics analysis of Hydra A. We note that there are no 
publicly available data of stellar kinematics for A478. Since 
the Hydra A cluster is at a very low redshift, strong-lensing 
information is not required for our joint analysis. 

We here summarize the dynamical properties of the 
BCG (3C218). It is well known that the AGN in the 
BCG ejects a pair of jets, sho wing interaction betwee n 
the jets an d surr ounding ICM llMcNamara et ahl l2000l ). 
iFuiita et al.l (l2013l) have discovered a dust lane in the cen¬ 
tral region of the BCG (4 kpc x 0.8 kpc). The major 
axis is along th e east-west di r ection , nearly perpendicular 
to the jet axis. iHamer et al.l (I20l4 ) found that the dust 
lane is a rotating disk by measuring several line emis- 
sions. The maximum rotating velocity reache s ~ 400 kms“^. 
iMelnick. Gopal-Krishna fc Terlevich I^1997^ also found a 


fast rotating disk using the Ho? line. Ekers fc SimkinI (Il983l ) 
measured stellar velocities and pointed out that the BCG 
is composed of a rotating disk of stars embedded in non¬ 
rotating stellar components. The rotational velocity inside 
4" (^ 4kpc) is found to be ~ 50kms“^ dLoubser et al.1 
1200^. The stell ar vel ocity dispersion is ~ 200-300 km s“^ 
( Heckman et al.l 1 19851 : iLoubser et al.ll2008l ). Therefore, the 
central part of the BCG ( ^ A") comprises cold rotating 
disk, rotating stellar, and non-rotating stellar components, 
exhibiting complex dynamics. On the other hand, the non¬ 
rotating stellar components dominate outside the rotating 
disk ( ^ 4”). 

We use the PIEMD model (Equation O, assuming 
spherically symmetry, isotropic orbits, and mass tracing 
light. We determine the core and scale radii by fitting the 
light distribution (Figure [5]) and then estimate the nor- 
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malization of the PIEMD and the NFW parameters by 
a joint fit to weak-lensing and stellar velocity dispersion 
measurements. The combined 'x function is defined by 
=Xg + Xdyn, where 

(9) 

^CTv* 


with and the velocity dispersion and its uncer¬ 
tainty, respectively. The velocity dispersion along the line- 
of-sight, (71 .0.3, is computed by the spherical Jeans equation 
(eq. 42 of ICappellari|[2008l ) 


S*(7l^o.s(^) — 2G 



j7*(r)M(r)F(r) 
-- dr 


( 10 ) 


(see also eq 12 of iNewman et al.l f201,‘ll 'l. where M{r) — 
MpiEMD,*(r) -I- MNFw(r) is the total mass, /3 = 1 — (jIIoI 
presents the velocity anisotropy, and p* and E* are the 
spherical and surface density profiles of the stellar compo¬ 
nents, respectively; here we use the best-fit light distribu¬ 
tion of the PIEMD model to describe i/* and E*. The kernel 
E'(r) is (r^ — for the isotropic-o r bit case, S = 0 (see 

ICappelIrf2008l : lElfasdottir et al.l2007l 'l. iLemze et al.l (l2012l 'l 
studied the radial dependence of /3(r) for simulated cluster¬ 
sized dark-matter halos, finding /3 ~ 0 — 0.2 at O.OSrvir. 
Since the line-of-sight integral is weighted by the stellar den¬ 
sity and thus dominated by the innermost region, the results 
are insensitive to /3(r) around the BCG. Hence, we can rea¬ 
sonably assume a constant /3. We fix /3 = 0 throughout the 
paper. We also checked that these results are not sensitive 
to the choice of j3 and do not change significantly when fit¬ 
ting with P = 0.2. We use the stellar velocity dispersion at 
r > 4" where the contribution from the non-rotating com¬ 
ponent is negligible. We summarize the lensing data, pro¬ 
jected velocity dispersion and the light profile for the BCG 
in Table [XT] The best-fit radial profiles of weak-lensing and 
projected velocity dispersion are shown in the left panel of 
Figure [T] The best-fit profile for the total lensing signal at 
r 100 kpc is systematically lower than the measurements, 
although consistent within errors. 

The spherical mass profiles M{< r) in the central re¬ 
gion (r < 300 kpc) are shown in the right panel of Figure 
[3 The stellar mass component (green dotted line) domi¬ 
nates over the NFW halo component (blue dotted line) at 
r 20 kpc. The NFW mass component from the joint fit is 
consistent with the best-fit NFW profile from the g+ fitting 
using the data at r > 150kpc (yellow region). The total 
mass profile (black solid line) agrees with the weak-lensing- 
only results using the full-range data (r > 30 kpc; cyan re¬ 
gion) because these best-fit tangential profiles are similar 
at r 100kpc (Figures |4] and (3). The NFW-|-point total 
mass profile from tangential-shear fitting (magenta region) 
is in good agreement at r 150 kpc with that from the joint 
analysis, but overestimates the joint results at r 50 kpc. 
We shall discuss this point in the last paragraph of this 
section as well as in Section 031 We note that the best-fit 
BCG mass profile is extrapolated beyond the maximum ra¬ 
dius for the light-profile measurements (~ 30 kpc) assuming 
a constant mass-to-light ratio. If the mass profile for the 
BCG is sharply truncated at this radius, the BCG mass is 
Mbcg = l-Oto'.4 X IO^^Mq. 

For comparison, we overplot the dynamical masses es¬ 


timated from the c old gas disk component (blue diamonds; 
iHamer et alj 20 131 and hydro static masses from X-ray ob¬ 
servations I David et al ] l200ll. red circles). A detailed com¬ 
parison of weak-lensing and X-ray mass mea surements at 
mass o v erdensities A < 2500 are presented in lOkabe et al.l 
l|2014bl ~l. iHamer et al.l l|2014l ') assumed a rotating exponen¬ 
tial disk for the central cold disk and obtained dynamical 
mass estimates within 3 .2 kpc and 2. 2 kpc u sing the Ha and 
Paa lines, respectively. IPavid et al.l (120011 1 derived hydro¬ 
static mass estimates using Chandra data. 

Overall, our results agree well with independent mea¬ 
surements from dynamical masses estimated from the cold 
gas disk and X-ray hydrostatic mass estimates, over a wide 
range of radii from a few kpc to 300 kpc, irrespective of 
different physical conditions and nature of cold gas parti¬ 
cles, hot diffuse plasmas, and stars used as mass tracers. 
Although the effects of non-thermal pressure (e.g., sound 
waves, subsonic gas motions, and cosmic rays) are expected 
to be significant as triggered by AGN activities, the X-ray 
and weak-lensing mass measurements in the central region 
are in excellent agreement. Since turbulent and/or rotational 
gas motions in the ICM can be directly proved by the SXS 
onboard Astro-H, a complementary combination of kinetic 
information about the ICM, cold gas, and stars and weak- 
lensing data will be able to provide an improved understand¬ 
ing of the physical state of the cluster center. 

We also estimate the stellar mass of the BCG using K- 
band luminositi es as described in Se ction [4T] An empirical 
scaling relation jArnouts et al.llioo^ l gives Mbcg,* ~ Slf x 
10 ^^Mq, within a photometric aperture radius of 17 kpc. 
The errors represent a 15% sca tter in the stellar mass-to- 
light ratio l|Arnouts et al.|[2007h . The inferred stellar mass 
is in remarkable agreement with the BCG mass determined 
by our joint analysis (Figure JIL Using the theoretic ally- 
predicted scaling relation jLonghetti fc Saraccd l2009l . Ta¬ 
ble 2), we find Mbcg,* 4 — 6 x 10 ^^Mq for the S alpeter 
IMF. On the other h and, the Chabrier (IChabrierl l2003l l 
and Kroupa dKrou ^ I 2 OOIII IMFs give a stellar mass of 
~ 3 X 10^^Mq, which is slightly lower than our best-fit esti¬ 
mate. 

The best-fit tangential shear pr ofile at r 100 kpc is 
some what discrepant with the data. iNewman et al.l (l2009l . 
l2013h took into account a calibration factor, mwL = 
P+, obs/5+,true, for then lensing model (i.e., p+, model ^ 
mwL X P+, model), possibly caused by systematic shear cal¬ 
ibration and/or redshift errors, to solve t he discrepancy 
betwe en their lensing and dyn amical data. Newm an et al.l 
ll2nn9h found mwL = 0.80, and INewman et al.l I 2013^ used 
a prior mwL = 0.89 ± 0.05, meaning that their shear mea¬ 
surements are underestimated by 00 10%-20%, relative to 
their dynamical measurements. 

We repeated htting including an additional correction 
factor mwL as a free parameter and obtained mwL = 
3.42 ± 1.24, which is extremely higher than our shear cal- 
ibration uncert ainty (2 — 3%; lOkabe et al.l l2013l . l2014al : 
lOkabe fc Smithll2015f l. It is therefore unlikely that the dis¬ 
crepancy is primarily caused by shear calibration errors, sug¬ 
gesting alternative possibilities, such as the choice of mass 
models. We estimate lensing signals from member galaxies 
located in the gap between the BGG-kinematics and weak- 
lensing data regions (8kpc ^ r Py 30kpc), assuming a stellar 
mass-to-light ratio based on the Salpeter IMF. Including 
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their contributions, the best-fit lensing profile is enhanced 
by ~ 20%, which however is insufficient to account for the 
discrepancy. In the joint analysis, we assumed that the BCG 
stellar mass traces the light distribution, and did not ex¬ 
plicitly include BCG’s dark-matter contribution relative to 
the NFW form. If there is a significant and extended dark- 
matter contribution associated with the BCG, the apparent 
discrepancy may be explained. We shall examine this hy¬ 
pothesis in the next subsection. 


4.3 Additional Dark-Matter Component around 
the BCG in Hydra A 

For the Hydra A cluster, the total mass profile around the 
BCG has been tightly constrained by the joint analysis of 
weak-lensing and dynamical data. We find that the inferred 
stellar mass contribution is not sufficient to fully explain the 
observed lensing signal g+ oc r~^ at r ~ 100 kpc. This may 
imply that there is an additional dark-matter component 
associated with the BCG. Here, we adopt as a working hy¬ 
pothesis that the central dark-matter distribution pDM{r) is 
shallower than that of the stellar component p*(r), namely 
PDM(r) < p*(r) at r ^ 0 and pDM(r) > p*(r) at larger 
radii. We have first attempted to fit with the generalized 
NFW model p oc (r/rs)~^(l -f and the cored 

NFW model p oc (1 -I- r/rc)~^(l -I- r/rs)~^, both described 
by three parameters, to account for the hypothetical dark- 
matter contribution. We find both models fall short of the 
central tangential shear signal and fail to explain the dis¬ 
crepancy. 

Now we employ the PIEMD model to describe the ad¬ 
ditional mass component around the BCG. We express the 
total mass in Equation (HOI as M(r) = MpiEMD,*(r) -|- 
MpiEMD,DM(r) -I- MNFw(r), where MpiEMD,DM(r) is spec¬ 
ified by three parameters, namely the normalization, core 
and scale radii. We here assume the minimum core radius 
to be 0.4' in order to satisfy the requirement that the dark 
matter distribution is shallower than the stellar distribution, 
because there is no visible tracer of hypothetical dark-matter 
profile, that is, the current data is not sufficient to constrain 
the shape of the profile. 

Figure |8] shows the resulting best-fit tangential shear 
profile g+{0) (left) and the spherical mass profile M{< r) 
(right). We find that adding the extended central dark- 
matter component can fully account for the the dynamical 
and lensing data. The change on the total mass profile due to 
the additional dark-matter component is small (right panel 
of Figure [8]). The resulting total mass profile agrees better 
with those from independent measurements. In Figure [9] we 
show the spherical mass density profiles for the total, dark- 
matter, NFW and BCG stellar components in the central re¬ 
gion. The dark-matter density profile at 1 kpc ^ r ^ 20kpc 
is shallower than the stellar density profile, as imposed by 
the prior (the bottom panel in Figure 0. The slope for the 
total matter density is perturbed by the BGG stellar and 
the additional dark-matter components from the slope of 
the NFW main halo. The minimized x^i„(d.o.f.) = 6.1 (20) 
has been improved from 9.0 (23) relative to the NFW model. 
The F test gives a probability of 5 x 10~^, thus supporting 
the addition of the flat dark-matter component to explain 
better both lensing and kinematics data. 

Lensing systematics, such as the effects of halo/BGG 


6 [arcmin] 

10 '^ 10 ° 



Figure 9. Top: mass density profiles in the central region of 
Hydra A. The black-solid, green-dashed, blue-dotted and red- 
dotted lines show the total, BCG stellar, dark-matter, and NFW 
mass density profiles, respectively. Bottom: mass density slope, 
—4 In p/d In r. The best-fit dark-matter density profile is shallower 
than the stellar density profile at Ikpc /S r 20kpc, as imposed 
by the prior and favored by the data. 


triax i alitv fe.g. [Gorless fc Kind l2 007l: lOeuri fc B landfordI 


I2OO9I : iMeneghetti et ah 2010l : Umetsu et al. 2015ah . cannot 
also be ignored as one of the possible causes of the discrep¬ 
ancy between weak-lensing data and the best-fit model for 
the NFW and BCG stellar mass components (Section 14.211 . 
Since the additional dark-matter component would not be a 
unique solution to explain the discrepancy, it is also crucial 
to investigate other possibilities. Sophisticated model of the 
triaxial/rotational/non-axisymmetric halo model for stellar 
kinematics and weak-lensing data would be critically impor¬ 
tant for a deeper understanding of the connection between 
the BGGs and host dark matter halo. Further systematic 
studies of very nearby clusters are of vital importance for 
understanding of the mass distribution around BGGs. 


4.4 Dark-Matter Fraction around the BCG in 
Hydra A 

Isobaric cooling time scales for the two clusters are less than 
1 Gyr within the central 30 kpc (iDavid et al.ll200l] : ISun et al.l 
l2003h . Gas cooling thus causes condensation of baryons, 
leading to formation of stars and contraction of dark mat¬ 
ter in the central cluster region. One of the most plausi¬ 
ble scenarios for the dark-matter resp onse to the cooling o f 
baryons is adiabatic contraction (e.g. iGnedin et al.|[2004l f. 
A joint analysis of stellar kinematics ( ^ 8kpc) and weak 
lensing data ( ^ 30kpc) for the Hydra A cluster allows us to 
decompose the observed total projected mass profile into the 
stellar and dark-matter components around the BGG, which 
provides us with a good opportunity to test the validity of 
adiabatic contraction models in the central region. 
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Figure 7. Left: tan gential reduced-shear f or th e Hydra A clus t er and projected velocity dispersion for the central BCG (inset). Circles 
and triangles denote iHeckman et al. and lLoubser et al] ll2008lj for projected velocity dispersion, respectively. The red-solid lines 

represent the best-fit models for the total lensing signal and the projected velocity dispersion. The dotted-green and blue-dot-dashed 
lines show the lensing signals for the BCG and NFW mass components, respectively. Right: spherical mass profile within 300 kpc. The 
back-solid and dashed-lines represent the total mass profile and its Icr uncertainty, respectively. The green-dotted and blue-dot-dashed 
lines represent the BCG stellar and NFW halo components, respectively. The magenta, cyan, and yellow shaded areas (from the top to the 
bottom) show the Icr error regions for the NFW-hpoint, NFW (r > 30kpc), and NFW (r > 150 kpc) models from the weak-lensing-only 
fits, resp ectively. The blue d iamonds, red circles, and green squares show dynamical m ass estimates from radio-line measurements of 
cold gas llHamer et alJ[2014^ . hydrostatic mass estimates from X-ray emitting hot gas llDavid et al.|[200lh . and stellar mass estimates 
from A'-band photometry assuming the Salpeter IMF, respectively. The top arrows indicate the radial ranges for the dynamical data 
and weak-lensing data. 
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Figure 8. Left: results for Hydra A including an additional extended dark-matter component associated with the BGG. The red-solid 
lines show the best-fit models for the total lensing signal and the projected velocity dispersion. The green-dotted, magenta-dashed, and 
blue-dot-dashed lines are the lensing signals for the stellar component, the dark-matter component associated with the BCG, and the 
NFW halo component, respectively. Right: spherical mass profile within 300 kpc. The back-solid, green-dotted, and blue-dotted-dashed 
lines show the results for the total mass, stellar mass, and total dark-matter distributions, respectively. For comparison, the total mass 
profile in Figure [7| is shown by the magenta-dashed line. The blue diamonds, red circles, and green squares are the same as those in 
Figure [8] 
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Here we estimate dark-matter mass fractions around 
the BCG, given by /dm = Mdm/(Mdm + Mt + Mgas), 
where Mdm is the dark-matte r mass , M« is t he stellar 
mass, and Mgas is the gas mass llPavid et al.ll200lh . We em¬ 
ploy the NFW (Sections 14.11 and 14.21) and composite dark- 
matter (NFW-fPIEMD; Section 14.31) models to calculate 
Mdm(< r). The stellar mass is composed of the BCG and 
cluster elliptical galaxies. The BCG stellar mass profile is 
calculated using the best-fit parameters given in Section r4.2l 
The stellar masses for other elliptical member galaxies are 
estimated by their Jf-band luminosities. Here, we assume 
the Salpeter IMF (Section 14.21) because the stellar mass es¬ 
timated by the Salpeter IMF agrees with that estimated by 
the stellar kinematics. 

Figure [TO] shows the resulting dark-matter fraction 
profiles for the Hydra A cluster. The cumulative dark- 
matter fractions progressively increase with increasing the 
cluster radius. The slope is slightly steeper beyond the 
effective radii of BCGs, ~ 20kpc, and flattened at 
r ~ 100 k pc. Similar results have been reported by earlier 
work ( e.g.lNag ino Sz Matsushitj 20091 : lOguri. Rusu fc Falcd 
I 2 OI 4 I '). iNagino fc Matsushita! (l2009l ~) found that the inte¬ 
grated mass-to-light ratio profiles for early-type galaxies 
based on the X-ray and photometry analysis increases be¬ 
yond effective radii of galaxies. lOguri. Rusu fc Falcd 
(l2014h conducted a joint analysis of strong lensing and stel¬ 
lar photometry for a large ensemble of elliptical galaxies, 
finding that the slope of the dark-matter faction becomes 
stepper at r > re. 

Now we examine the dark-matter fractio ns using mod¬ 
els w ith and without adiabatic c ontraction (iGnedin et all 
l2004h . We use the Jaffe model lljaffel Il983l') as a stellar 
mass distribution, where the density profile follows p(r) oc 
r~^{r + s)~^. In the limit of r ^ a for the PIEMD 
model (Equation (O), the PIEMD density profile resem¬ 
bles the Jaffe profile For the dark-matter component, we 
only consider the NFW model. The dark-matter fractions 
dashed lines) pred icted by adiabatic contraction models 
Gnedin et al.l 120041 ') overestimate our measurements. On 
the other hand, models without adiabatic contraction agree 
better wit h the dark-mat t er fra ctions , consistent with the 
result s of I Newman et al.l l|201,3h and lOguri. Rusu fc Falcd 
l|2014h . This suggests that other physical processes are criti¬ 
cally important to suppress the modification of dark-matter 
profile. 


5 DISCUSSION 

We fi nd that the adiabatic contraction model (iGnedin et al.l 
l2004h overestimates the dark-matter fraction profiles com¬ 
pared to our measurements (Section 14.41) . One of possible 
scenarios is that the gas outflows driven by AGN activities 
suppress condensa tion of the dark-matter distribution in the 
central region (e .g. iTevssier et aI]l201ll:r Martizzi et al.l[201^ : 

iRagone-Figueroa. Granato fc Abadill2012[) . In fact, on-going 

jet ac tivities of central AGNs have been reported in two clus- 
ters ([David et al.l 1200 ll : ISun et all l2003l : IPiehl et al.l I2OO8I . 
e.g.). The projected distances between the BGG and X- 
ray bubbles triggered by the observed activities are r ~ 
25 kpc for Hydra A and r ~ 9 kpc for A 4 78, respec- 
tively (e.g. I David et al.l I2001I : ISun et akl I2OO3I : IWise et ahl 



Figure 10. Dark-matter fractions as a function of spherical clus¬ 
ter radius. The red circles and blue triangles show the results for 
Hydra A using the best-fit NFW (Section l4.2ll and NFW+PIEMD 
(Section l4.3l l profiles, respectively. The solid and dashed lines are 
the dark-matter fractions predicted by models without and with 
adiabatic contraction, respectively iGnedin et ahlEoOlh . The two 
arrows mark the radial range which covers each data. 


I 2 OO 7 I : IPiehl et al.l l2008h . Three-dimentional characteristic 
scale radii obtained by equating the interior stellar and 
NFW masses are req ~ 25 kpc for Hydra A; Teq ~ 
22 kpc and ~ 15 kpc for A478, assuming the Salpeter 
and Chabrier IMFs, respectively. This suggests that the 
current AGN jet activities can directly affect the baryons 
over the region where stellar mass dominates. Although 
there is an uncertainty between the three-dimentional ra¬ 
dius, req, and the projected distances between the BCG and 
X-ray bubbles, the results hold even if the jet/bubble di¬ 
rection is ^ 40° from the plane of the sky. Another alter¬ 
native scenario assumes mi nor, dry (dissipationless) merg¬ 
ers of clumpy structures (e.g.lEl-Zant. Shlosman fc HoffmanI 
I 2 OOII : iLackner fc Ostrikerll201o[l . which can flatten the cen- 
tral dark-matter density profile (Section 14.31) and be closely 
correlated with the BCG growth. 

Further systematic studies of central dark-matter mass 
fractions and density slopes, using a large sample of AGN 
and non-AGN clusters, are required to distinguish between 
the effects of gas outflows and dissipationless mergers, for 
deep understanding of physical processes governing the BCG 
growth. Systematic joint studies combining weak lensing and 
stellar kinematics data will also be useful to examine corre¬ 
lations between stellar population p roperties and IMF con - 
straints inferred from joint analyses (iMcDermid et alll2014l l. 


6 SUMMARY 

We have presented a weak-lensing study of the nearby cool- 
core clusters Hydra A and A478. To minimize dilution of 
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the lensing signal due to contamination by cluster mem¬ 
bers, we carefully select populations of background source 
galaxies in the color-magnitude plane. Compared to weak- 
lensing studies of c lusters at intermediate redshifts 2 ~ 0.2 
llOkabe et al.ll2013 '). the level of background contamination 
by cluster members is significantly reduced. This enables us 
to obtain a larger number of background galaxies for weak- 
lensing measurements while achieving a low level of con¬ 
tamination. This is one of the advantages of weak-lensing 
analysis of very nearby clusters. 

The S/N ratios for detection of the tangential shear sig¬ 
nals are 7.0 and 8.1 for Hydra A and A478, respectively, after 
accounting for the contribution from projected uncorrelated 
LSS. The resulting S/N ratios are comparable to those of 
clusters at 2 ~ 0.2, thanks to the increased number of back¬ 
ground galaxies behind the clusters. 

We find that the tangential shear profile for A478 is well 
fitted with a single NFW model. For Hydra A, the tangen¬ 
tial shear signal in the central region is well described by 
a two-component model including the central BCG as an 
unresolved point mass. 

We find that the choice of IMFs can introduce a large 
uncertainty (factor of ~ 2) in the BCG stellar mass esti¬ 
mates, which makes it difficult to decompose the observed 
total projected mass profile into the stellar and dark-matter 
component s CSection 14.111 . On the other han d, as demon¬ 
strated by iNewman et al.l (l2009l . I 2 OIII . l2013l'l , the internal 
stellar kinematics of BCGs enables us to precisely measure 
the mass profiles of two components independent of the IMF 
uncertainty fSection l4.2l) . For Hydra A, we find that the cen¬ 
tral mass profile (< 300 kpc) determined from weak lensing 
is in excellent agreement with those from independent mea¬ 
surements, including dynamical masses estimated from the 
cold gas disk component. X-ray hydrostatic total mass esti¬ 
mates, and central stellar mass estimates. For the BGG, the 
data prefer the Salpeter IMF to the Chabrier and Kroupa 
IMFs. An additional flat dark-matter component around 
the BCG accounts simultaneously for the weak-lensing and 
stellar-kinematics data for Hydra A. 

Dark-matter fractions around the BGGs for the Hy¬ 
dra A cluster are found to be smaller than those predicted 
by adiabatic contraction models, and to agree well with 
model predictions without contraction. This implies that 
other baryonic processes, such as the AGN feedback, dis¬ 
sipationless mergers, or the combination of the two, could 
play important roles in shaping the central cluster mass pro¬ 
file. A precise joint measurement of the central cluster mass 
profile provides us with complementary information about 
the dynamics of hot intracluster gas, which will be directly 
probed by Astro-H. Stellar kinematics data for the BCG of 
A478 is also important for future works. 
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APPENDIX A: MULTIWAVELENGTH DATA 
OF THE HYDRA A CLUSTER 

We here summarize the observing data and their measure¬ 
ment errors for the tangential shear profile, the projected 
velocity dispersion and the light profile for Hydra A. 

This paper has been typeset from a LjijX/ DTf^X file prepared 
by the author. 
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Table Al. Tangential shear {g+), projected velocity dispersion 
(cTt;,*) and the core (a) and scale radii (s) of the PIEMD for the 
Hydra A. Note that the measurement error for the tangential 
shear { C = C*stat + C*t., S s) is weakly correlated. ^ iHeckman et af] 
lll985lj .^ lLoubser et alj ll2008lj . 


e 

9+ 

Cl/2 

(arcmin) 



5.65922 X 10-1 

1.75256 X 10-1 

1.12041 X 10-1 

7.13965 X 10-1 

1.28711 X 10-1 

7.97313 X 10-2 

9.03514 X 10-1 

7.36090 X 10-2 

6.06401 X 10-2 

1.12847 

7.98743 X 10-2 

4.91149 X 10-2 

1.43289 

5.48124 X 10-2 

4.17630 X 10-2 

1.79821 

5.15705 X 10-2 

2.83264 X 10-2 

2.27188 

4.08415 X 10-2 

2.18615 X 10-2 

2.83350 

2.58124 X 10-2 

1.76992 X 10-2 

3.57655 

2.13949 X 10-2 

1.37130 X 10-2 

4.51228 

7.67107 X 10-3 

1.09288 X 10-2 

5.69914 

1.69769 X 10-2 

8.85021 X 10-3 

7.14476 

7.01793 X 10-3 

7.00017 X 10-3 

8.97936 

1.40250 X 10-2 

5.61649 X 10-3 

1.12947 X 10 

1.04568 X 10-2 

4.69525 X 10-3 

1.42627 X 10 

9.69863 X 10-3 

3.89583 X 10-3 

1.78611 X 10 

7.65808 X 10-3 

3.38880 X 10-3 

2.25777 X 10 

7.34974 X 10-3 

2.91867 X 10-3 

2.84133 X 10 

4.07282 X 10-3 

2.49727 X 10-3 

3.56222 X 10 

1.74820 X 10-3 

2.24314 X 10-3 

4.41888 X 10 

2.57557 X 10-3 

2.20707 X 10-3 

e 



(arcmin) 

(kms-^) 

(kms-^) 

6.83333 X 10-2t 

3.29000 X 102 

5.70000 X 10 

7.00000 X 10-2t 

3.83000 X 102 

6.10000 X 10 

8.03333 X 10-2t 

3.29000 X 102 

2.60000 X 10 

9.50000 X 10-2t 

2.54000 X 102 

1.18000 X 102 

1.01667 X 10-it 

2.85000 X 102 

8.20000 X 10 

1.30000 X 10-it 

3.39000 X 102 

1.13000 X 102 

a 

s 


(arcmin) 

(arcmin) 


1.751265 X 10-2 

4.68178 X 10-1 



















